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ABSTRACT 

We have studied non-stationary and non-axisymmetric perturbations of a magneto- 
hydrodynamic accretion onto a rotating (Kerr) black hole. Assuming that the magnetic 
field dominates the plasma accretion, we find that the accretion suffers a large radial 
acceleration resulting from the Lorentz force, and becomes highly variable compared 
with the electromagnetic field there. In fact, we further find an interesting perturbed 
structure of the plasma velocity with a large peak in some narrow region located slightly 
inside of the fast-magnetosonic surface. This is due to the concentrated propagation 
of the fluid disturbances in the form of fast-magnetosonic waves along the separatrix 
surface. If the fast-magnetosonic speed is smaller in the polar regions than in the equa- 
torial regions, the critical surface has a prolate shape for radial poloidal field lines. 
In this case, only the waves that propagate towards the equator can escape from the 
super-fast-magnetosonic region and collimate polewards as they propagate outwards in 
the sub-fast-magnetosonic regions. We further discuss the capabilities of such collimated 
waves in accelerating particles due to cyclotron resonance in an electron-positron plasma. 

Subject headings: acceleration of particles — accretion — galaxies: active — MHD - 
relativity 
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1. Introduction 

It is commonly accepted that extragalactic jets are 
intimately linked with accretion process onto super- 
massive black holes residing in the central regions of 
active galactic neuclei. These jets are initially rela- 
tivistic, as indicated by superluminal proper motions 
of radio emitting knots (e.g., Wehrle et al. 1992) and 
by high-energy, rapidly variable 7-ray emissions (e.g., 
Montigny et al. 1995). Moreover, Hobble Space Tele- 
scope studies of the base of M87 jet reveal a rotating 
gas disk apparently lying normal to the jet direction 
(Ford et al. 1994; Harms et al. 1994). Despite intense 
study, the underlying formation mechanism is still un- 
certain. Nevertheless, hydrodynamic and magnetohy- 
drodynamic processes associated with the accretion 
disk seem to be a promising candidate mechanism. 

Possible flows of the energy conversion from the ac- 
cretion to a small fraction of gas in jets have been sug- 
gested by Shakura and Sunyaev (1973) in the context 
of a thick supercritical accretion disk which exhib- 
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poles. If the radiation-supported rotating gas adopts 
a hydrostatic toroidal configuration, then a pair of 
funnels are defined which could be responsible for the 
production of jets along the rotational axis (Lynden- 
Bell 1978). Blandford and Payne (1982) considered 
a magnetized disk and showed that a gas leaves the 
disk in a centrifugally driven wind provided that the 
magnetic field makes an angle of less than 60° with 
the radius vector at the disk. Furthermore, the mag- 
netohydrodynamic (MHD) disturbances produced at 
a galactic nucleus with a compact nuclear disk would 
be strongly collimated polewards resulting in jets, if 
the Alfven velocity in the disk is much higher than its 
surroundings (Sofue 1980). Such a collimation would 
lead to an increase in the wave amplitude, resulting 
in shock waves which are conjectured to develop into 
a strongly compressed region of the magnetic field. 
In this region, high-energy particles are likely to be 
accelerated in the perpendicular direction to the equa- 
torial nuclear disk. 

The purpose here is to explore the issue whether 
MHD waves can convey some portions of accretion 
energy to the polar regions in the vicinity of a black 
hole. Causality requires that the MHD inflows should 
pass through the fast-magnetosonic point and become 
super-fast-magnetosonic at the horizon. The investi- 
gation of the so-called critical condition that the in- 
flow should pass through this point smoothly offers, 



in fact, the key to an understanding of MHD interac- 
tions in a black hole magnetosphere. 

The MHD interactions are expected to work most 
effectively in the magnetically dominated limit in 
which the rest-mass energy density of particles is neg- 
ligible compared with the magnetic energy density. In 
this limit, the fast-magnetosonic point is located very 
close to the horizon (e.g., Phinney 1983); as a re- 
sult, a general relativistic treatment is required. An- 
alyzing the critical condition in a stationary and ax- 
isymmetric magnetically dominated black hole mag- 
netosphere, Hirotani et al. (1992, hereafter Paper I) 
showed that roughly 10% of the rest-mass energy and 
a significant fraction of the initial angular momentum 
are transported from the fluid to the magnetic field 
during the infall. Furthermore, if a small-amplitude 
perturbation is introduced into the magnetosphere, a 
lot of perturbation energy is deposited from the mag- 
netic field to the fluid near to the fast-surface in the 
short- wavelength limit; accordingly, the plasma ac- 



crction becomes highly variable there (Hirotani et al 



1993 , hereafter Paper II) . Subsequently, Hirotani and 



Tomimatsu (1994, Paper III) investigated the spatial 
structure of the disturbance in a Schwarzschild met- 
ric, assuming that the characteristic scale of the ra- 
dial variations of perturbed quantities is comparable 
to that of unperturbed quantities instead of adopting 
the short- wavelength limit. They revealed that the 
magnetically dominated accretion becomes most vari- 
able at the fast-magnetosonic separatrix surface and 
that the large-amplitude fluid's disturbance can es- 
cape into the sub-fast-magnetosonic regions by prop- 
agating meridionally almost along the separatrix. In 
this paper, we extend the analysis performed in Paper 
III to a Kerr metric, further examine the propagation 
of the escaped fast waves in the sub-fast regions, and 
discuss the possibilities of particle acceleration due to 
a collimation of such waves. 

The outline of this paper is as follows. In §2, 
we formulate non-stationary non-axisymmetric per- 
turbations of the MHD accretion and derive the wave 
equation describing the perturbation. Solving the 
wave equations, we show in §3 that the fluid becomes 
most variable slightly inside of the fast-surface, which 
is consistent with the results obtained in Paper II and 
III. We further demonstrate that the disturbances can 
escape into the sub-fast regions in the form of fast 
waves by propagating towards the equator provided 
that the fast-magnetosonic speed, C/fm is slower in the 
polar regions than in the equatorial regions. In §4, the 
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escaped fast waves will be shown to collimate towards 
the rotational axis under such a distribution of Ufm- 
We finally discuss the capabilities of such collimated 
waves in accelerating of particles due to nonlinear in- 
teractions between the waves and electron-positron 
plasmas in §5. 

2. Magnetically Dominated Accretion 

We will begin by considering basic equations de- 
scribing a magnetosphere around a rotating black 
hole. Since the self-gravity of the electromagnetic 
field and plasma around a black hole is very weak, 
the background geometry of the magnetosphere is de- 
scribed by the Kerr metric, 



ds 2 = 



A - a 2 sin 
£ 

Asm 2 9 , 



-dt 2 



4Mar sin 2 9 



-dr 2 



E 



dtd(j) 



(1) 



I A 



where A = r 2 — 2Mr + a 2 , £ = r 2 + a 2 c 
(r 2 + a 2 ) 2 - Aa 2 sin 2 9 and a = J/M; M is the mass 
of a hole. Throughout this paper we use geometrized 
units such that c = G = 1 . 

Under ideal MHD conditions, since the electric 
field vanishes in the fluid rest frame, we have F^U^ = 
0, where is the electromagnetic field tensor sat- 
isfying Maxwell equations, F[^ v<p ] — and C/ M is the 
fluid four velocity. The motion of the fluid in the cold 
limit is governed by the following equations of motion: 



2.1. Unperturbed Magnetosphere 

From an analysis of the stationary and axisymmet- 
ric ideal MHD equations, it is known that there exist 
four integration constants that are conserved along 
each flow line (e.g. 



Bekenstein and Oron (1978): Ca 



mcnzind 1986a, b). These conserved quantities are the 
angular velocity of a magnetic field line (Qf), the par- 
ticle flux per magnetic flux tube (77) , the total energy 
(E) and the total angular momentum (L) per particle. 
They are defined as follows: 



Op = 



F t e 



rt/i 



(3) 



-gnU r 



-gnU 6 



F e $ 



and 



E = fxU t - ^B 4 

4:717] 



(4) 



(5) 



(G) 



where the toroidal magnetic field (B^) is defined by 
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B& — — 



-Fri 



-9 



(7) 



[mU»U" + ^-{F»PF p v + X -g» v F a!i F af} ) 



= 0, 



(2) 



where the semicolon denotes a covariant derivative 
and fi the rest-mass of a particle. For electron-proton 
plasmas /z refers to a rest-mass of a proton, whereas 
for electron-positron plasmas /x refers to that of an 
electron (or a positron). The proper number density 
(n) obeys the continuity equation, (nU^)-^ = 0. We 
adopt these basic equations for a description of sta- 
tionary and axisymmetric black-hole magnetospheres 
in §2.1 and for an analysis of perturbed state in §2.2 
and afterwards. 



2 = Asm 2 9, 



1 — g = £ sin ( 



(8) 



When < ftp < 51h holds, thereby hole's rotational 
energy and angular momentum are extracted mag- 
netohydrodynamically (Blandford and Znajek 1977; 
Phinney 1983), both E and L become negative. A 
poloidal flow line is identical with a poloidal magnetic 
field line, and is given by ^>(r, 9) = constant, where 
^ is the 0-component of the unperturbed electromag- 
netic vector potential. The conserved quantities are 
functions of 'J alone. 

We next describe a stationary plasma accretion in 
a black-hole magnetosphere. In a black-hole magne- 
tosphere there are two light surfaces. One is called 
the outer light surface, which is formed by centrifu- 
gal force in the same manner as in pulsar models. 
The other is called the inner light surface, which is 
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formed by the gravity of the hole. In a region be- 
tween the horizon and the inner light surface, the 
plasma must stream inwards, while in a region be- 
yond the outer light surface it must stream outwards. 
A plasma source in which both inflows and outflows 
start with a low poloidal velocity is located between 
these two light surfaces (Nitta et al. 1991); the injec- 
tion region (r = n) of the accretion may be a pair- 
creation zone above the disk (Beskin et al. 1992; Hi- 
rotani and Okamoto 1997) or the disk surface whose 
inner edge corresponds to the innermost stable cir- 
cular orbit in a Kerr geometry. Along the magnetic 
field lines the plasma inflows pass through the Alfvcn 
point, the light surface, and the fast magnetosonic 
point (r = tf), successively; they finally reach the 
event horizon (r = ru). From now on we use the sub- 
scripts I, F and H to indicate that the quantities are 
to be evaluated at r — n, r = rp and r — rn, respec- 
tively. In a magnetically dominated magnetosphere, 
the fast-magnetosonic point is located very close to 
the horizon (for explicit expression of rp and rjj see 
Paper I). 

2.2. Perturbed Magnetosphere 

We next consider a small-amplitude non-axisymmetric 
perturbation superposed on the unperturbed state 
discussed in the last subsection. In the perturbed 
state all perturbation equations are solved self-consistently, 
including the trans-field equation. We wish to ex- 
amine the behavior of fluid quantities (energy, angu- 
lar momentum, and poloidal velocity) in response to 
small variations in the magnetic field at various points 
in the magnetosphere, especially at the fast-point. 

Let the actual poloidal component of the electric 
and magnetic fields, as a result of the disturbance, 
be Ea + e-A and Ba + ^a {A = r, 9) , respectively; 
the small letters (e r , eg, b r , and be) are the Eulerian 
perturbations of the corresponding quantities. Let 
and &0 denote perturbations of F tt f, and —^—gF r e, 
respectively. Furthermore, we introduce u r and u 9 
such that the actual component of the poloidal fluid 
velocity field, as a result of the disturbance, becomes 
U A + u A . Let ut and — u^, be the Eulerian perturba- 
tions of fluid energy and angular momentum per unit 
mass, respectively. 

Making use of (4), we can simplify the t, <f> compo- 
nents of the frozen- in conditions as follows: 



U*e^ + V r yf~gb B - ^ g u e = 0. 



(10) 



The (^-component of the frozen-in condition in a little 
bit complicated. To take account of cancelations in 
the A-expansion, we combine U^u^ = with the 9- 
component of the frozen-in condition. We then obtain 



r (Gie e + G 2 V^gb r ) + g rr geeU r 



= + ~7T U 4 
9 M G\ 



where 



-9 B <p[{t/ ^)U<i>- (gtg + ff00^F)] t 
(Gi+G 2 f! F )Gi 



Gi = g^Ut — gt<f>U$, 
G2 = —gttUj, + gt<t>Ut- 



0(11) 



(12) 
(13) 



Since we are concerned with magnetohydrodynamical 
interactions near to the fast-surface located very close 
to the horizon, we evaluate unperturbed quantities 
appearing in the perturbed equations in the vicinity of 
the horizon (A <C A/ 2 ). Then equation ( pi] ) reduces 
to 



2Afr H SH sin 2 6> 



A 

jjr 
H 



r \2 



2Mn 



sin 9 



-uf + (n H - r» F ) 5 H 



u r = 0, (14) 



where 



e 
h 



E 



n F L, 



(15) 
(16) 



fin is the hole's rotational angular frequency and de- 
fined by 



(17) 



Here use has been made of the fact that unperturbed 
fluid velocity t/g is related with fluid energy t/ t H and 
angular momentum — 1/¥ as (Paper I) 



U r e r + n F ^ e u e + £/% = 0, 



(9) 



05 - -^-(U? + SlnUf) + 0(A/M 2 ). (18) 
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UjJ and hence £/ t H and —U^, are of order of unity. 
The explicit expressions of these quantities are given 
in Paper I. 

If we assume that the characteristic scales of merid- 
ional variations in the perturbed state are much 
shorter than those in the unperturbed state, the per- 
turbation equations reduce significantly. For non- 
axisymmctric perturbations, all perturbed quanti- 
ties may therefore be assumed to be proportional to 
exp(iujt — ikgO — im(f>), where \kg\ 1. Under this 
approximation, we obtain from the homogeneous part 
of Maxwell equations 



imb^ + A sin 2 9 [ - //.-,;// 
2Mr H de 



2Mr H db r 
A dr* 



0, (19) 



A dr* 



+ ime r + itu^/^gb = 0, (20) 



ikge r + 



2Mr H de e 
A dr.. 



Asin( 



0, 



(21) 



Substituting equation (|2^) into equation ( pl| ) , we have 



dee 
dr* 



lu — mf2p 
. S H b 



de$ 
dr* 



-ilu- 



2Mr H sin 9 

In the same manner, equations ( |l9| ) and ( p^ ) give 



(25) 



-dF_ 



to ■ 



-im 



de^ 
dr* 



2Mr H sin0 



(26) 



Equations (|23|)-(|2q) express poloidal components of 
perturbed electromagnetic fields in terms of their 
toroidal components. 

Let us next consider the equations of motion. First, 
the definition of proper time U^u^ = yields near the 
horizon with the aid of equation (fL8h , 



where r* is the tortoise coordinate defined by 



dr* 
dr 



(22) 



It is convenient to introduce this coordinate when 
we describe waves near to the horizon, because the 
interval (rjj,oo) in the r-coordinate is stretched to 
(—00,00) in r*. We assume that the characteristic 
scale of the radial variations is comparable to that 
of the gravitational field, that is, \df/dr*\ w \f/M\, 
where / denotes some perturbed quantity. In this 
paper we do not take the short-wavelength limit 
(\df/dr*\ ^> \f/M\) which was adopted in Paper II. 

Eliminating u e from equations @ and (|Tb"|), we 
obtain a relation between e r , b e , and e^. Combining 
this relation with equation (po|), we can express e r 
and b e in terms of as follows: 



(to — mOfjer 



2Mn 



de<j. 
dr~ 



(lu — mf2p)-\/ —gb 6 



2Mn 



dej 
dr* 



im(Qn — £If)£ 



(23) 



(24) 



u r = _^H (ut + n J + 0(A/M 2j_ (27) 

Secondly, the ^-component of the equation of motion 
gives 



(imU»V»U t )V + fW = 0; 
in the vicinity of the horizon we have 



A = —i(to — mCln) + -; — , 
dr* 



(28) 

(29) 
(30) 



f(i) 



47TET 



A 



i(oj — mSinJ— ^; n - 



sin 9 



(31) 



The subscript (1) indicates that the quantity is eval- 
uated in the perturbed state. In general, the a- 
component of the Lorentz force is defined by 



_ F a ^(v=9 F n 



47I\ 



(32) 
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Setting a = t, taking the linear order in the pertur- 
bation, and evaluating in the vicinity of the horizon, 
we obtain equation 

The perturbed fluid's density, ni, can be elimi- 
nated from equation (pq) with the aid of the conti- 
nuity equation 



(33) 



where 



H(u) 



du r .. ,2Mr H , „ . 
h i{u - mil s ) — (u t + ilnUtp) 

-ike—— — v e . 
2Mr H 



(34) 



Operating A to the both sides of equation (|33|), and 
combining with (p8|) to eliminate m, we obtain 



r H -M 
Mm 



-A/« = 



A ) ut - (d r U t )H(u) 
(35) 



In the same manner, from the ^-component of the 
equation of motion, we obtain 



A 



+A/i 1} - 0, 



(36) 



where 



/ i 1) = "fV /tliJ+ ° (A/M " ) - (3?) 
Finally, ^-component of the equation of motion be- 
comes 



p(i) 



(38) 



where 



Eliminating n\ with the aid of equation (|33|), tak- 
ing the leading orders in A-expansion and consider- 
ing relative amplitude between perturbed quantities 
together with the dispersion relation for the outgo- 
ing fast-magnetosonic mode (Paper II), we can self- 
consistently neglect fluid contributions in equation 
(pq) to obtain 



A/f x) = 0. (40) 
i3|) , (BJ) , we can rewrite equa- 



lising equations (jl 
tion into 



(2Mr H ) 2 fi H - «f / r H - M 



AEh w - 



mf2p 
2Mr H 



2Mr r 



= 



\ I Ae 

(41) 



We analyze a system which is formed by 10 equa- 
tions (equations @, @, g-g, and 05)) 

in 10 unknown functions (b r , b\e ri eg, b#, u t , u r , 
u 9 , and zt^). The perturbed fluid density would be 
calculated if we solved the perturbed continuity equa- 
tion. These 10 equations give some relations between 
two arbitrary perturbed quantities, and are further 
combined into a single differential equation. 

Substituting equation ( pi[ ) into equation (|l(]), we 
obtain 



i^DE^Ae* (42) 



We can use this equation to eliminate u and examine 
the relative amplitude between u and e^. 

One combination of equations (|35|) and (|3^) yields, 
with the aid of equation (07]) , 



E Sh 



sm f 



(43) 



2Mr H (2Mr K 



A(e r + n^yf^b 8 ) - 



2Mr H 

ike B &0 u - m£l F 
2Mr H sin6» _ * 2Mr H sinJ 



e (39) 



Here, use has been made of the fact that the unper- 
turbed fluid quantity 

P ( +fl F ^ = ^^M (44) 

is constant along each flow line. The other combina- 
tion of equations (p33) and (pq) yields 
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u t = -n F u* + 0(A/M 2 ) 



(45) 



Using equation ( p7|) and ( f45|) , we can eliminate 
in equation ( |l4|) to obtain 



(51) 



and recovering meridional derivatives by setting kg = 
idg, we can rewrite equation (E9I) into 





A 


'2Mr H , 
/i 


6* 




sin# 



(46) 



where the effective potential K is defined by (Taka- 
hashi et al. 1990) 



K = g t t + 2g t 4,flF + 94,4,^%- 



(47) 



At the injection point where U r vanishes, K takes a 
positive value, K\ = (e//i) 2 , whereas at the horizon 
we have K = Kn = g^A^n - ^h) 2 < 0. Moreover, 
equations (^5| ) and ( |2q ) give 



dh_ 



(f^H - ^F)h 

u — mfip 



Ae ri 



-i(uj — mflu) 



2Mr H sin9 



(48) 



So far, we have obtained four independent equa- 
tions (|l]), d||), HI), and @ for four unknowns u r , 
/i = eg — flay/— gb r , e0, and 6^. Combining these four 
equations, we finally obtain the wave equation 



(A-A F 



dr 2 
r H - M 



Mm 
i{uo — mfifj) 
A ^ 2 



(A + A F ) + 2i(u - mO H )A F 

m-M 



fee 2 \u r = 0, 



+ i(w — mfiu) 



,2Mr H> 

where Ap is given by (Paper I) 



d 

(A + A F ) 

(49) 



A _ 2Mr H S H sin 2 0Q F (n H - » F )(t/fi) 3 M 



Kj-K-a 



E 



Introducing a new non-dimensional radial coordinate 
as 



x 2 (x — Xf) 



- u 



2 + 2x(x + ictxf)-^ 



dx 

+ ia)(x 



■ xf) 



dx 
Z H x 2 



d 2 



4(r H - M) 2 89 2 







(52) 



where non-dimensional corotational frequency a is de- 
fined by 



a — (oj — mflu) 



Mm 



(53) 



Equation ( p2j ) becomes elliptic in the sub-fast region 
(x > xf), while it becomes hyperbolic in the super- 
fast region [x < xf)- 



From equations ( |27[ ) and (j45|), we can see that 
fluid's energy (u t ) and angular momentum (— Uj,) 
obey the same equation as (|52|). Therefore, equa- 
tion (|52|) describes the fluid's disturbance near to the 
horizon. In the slowly rotating limit (a — > 0), equa- 
tions (|5l|)-(|5^) reduces to equation (25) in Paper III, 
in which axisymmetric (m = 0) perturbations in a 
Schwarzshild metric were examined. 

3. Highly Variable Accretion 

To examine the spatial structure of u r (x, 9), let us 
rewrite (^2|) as 



U 1 ■ d \ r, S H^ 2 d 2 \ n 

\{- 1 - ta + x m) DFM + 4(m-M) 2 W\ = ° 

(54) 

where Dfm refers to the differential operator asso- 
ciated with outgoing-fast-magnetosonic mode and is 
defined by 



Dfu = x(x — Xf)- — h ia(x + Xf) + x (55) 
ox 

Let us examine the radial(a;) dependence of u r by ne- 
glecting the 0-derivative term as a first step. Under 
this assumption, the ingoing and the outgoing modes 
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in equation ( p4[ ) can be completely separated. Equa- 
tion Dfmu t = gives a solution corresponding to the 
outgoing mode, 



C v 



+2ia ' 



(56) 



where C\ is an integration constant. In order that 
the right-hand-side may not diverge at x = Xp(6), 
the real part of 1 + 2io~ should be non-positive; this 
indicates that outgoing radial waves must decay, be- 
cause no steady supply of perturbation energy across 
the fast-surface x = xp(9) is possible. In this paper, 
we postulate a steady excitation of perturbation in- 
duced by plasma injection from the equatorial disk 
or a pair production zone above the disk. It requires 
that a should be real. Then, to avoid divergence at 
x = xp, we must consider the ^-derivative term in 



(54) at least near to the fast-surface. 

Let us modify the solution (|5l) into the form 



What is most important here is that this maxi- 
mum amplitude surface coincides with the separatrix 
of characteristics. To see this, it is convenient to in- 
troduce a new radial coordinate £ which denotes a 
deviation from the maximum-amplitude surface, 

£ ee x - x F + Re(<5) (61) 
The characteristics of equation (b3) are expressed as 



dx 
d9 



(62) 



In the super-fast region, any waves must propagate 
inwards, dx < 0. Therefore, waves propagating to 
lower latitudes (d9 > 0) are indicated by the upper 
sign, while those to higher latitudes are indicated by 
the lower sign. Combining ©, ©, and ©, we 
obtain an equation expressing how a characteristic 
deviates from the maximum-amplitude surface, 



x - x F (6) + 8(6)] 



(57) 



where 5 should be a complex function of 9, so that we 
may obtain a regular solution u r for real a. Inserting 
(57) into ( ]54| ) and evaluating the equation in the limit 
| (x — xp + e)/xp |< 1, we obtain a non-linear first- 
order differential equation for 8, 



( dx-p 
= \~dT 



dS 



(58) 



As boundary conditions, we impose the following 
symmetry conditions, 



dRe(S) _ dlm(S) 

d9 ~ dO 
dRe(8) dlm{8) 



dO 



dB 







at 
at 



0, (59) 
I (60) 



As an example, numerical solutions satisfying these 
conditions for xp = 0.1 cos 2 9 + 0.01 (i.e., prolate 
shape of the fast-surface) are depicted in Figs. |]and||. 
In Fig. |, Re(8) = {dx F /d9) 2 = 0.01 sin 2 (2(9), which 
approximately satisfy equation when d 2 Xp/d9 2 w 
(i.e., at 9 w tt/4), is depicted by the dashed line for 
comparison. It should be noted that the maximum 
amplitude surfaces, x — x F — Re(<5), where u r has a 
sharp peak, is located slightly within the fast surface, 
x = xp. 



dd 



dx-p I d9 
\dx F /d9\ 



^RH8). (63) 



It follows that for a prolate shape of the fast-surface 
(dxp/d9 < 0), dt; has the same sign as £ for waves 
propagating to lower latitudes (dO > 0, upper sign). 
In other words, outside of the maximum amplitude 
surface, characteristics deviate from this surface out- 
wards, whereas inside of this surface, they deviate in- 
wards. Poleward propagating waves (d9 < 0, lower 
sign), on the other hand, always deviate inwards 
(d£ < 0) to be swallowed by the hole. That is, only 
waves propagating towards the equator can escape 
into to sub-fast regions, if the fast-surface is prolate 
(Fig. |^) . Thus we can regard the maximum amplitude 
surface as separatrix of characteristics. The same dis- 
cussion could be applied for a oblate shape of the 
fast-surface. 

Since fluid obtains most of its perturbation energy 
from the electromagnetic field at the maximum am- 
plitude surface, £ = 0, most of the fluid's disturbance 
propagates almost along the separatrix (i.e., the max- 
imum amplitude surface) and at last deviates inwards 
or outwards. In other words, meridional propagation 
is essential to examine the spatial structure of fluid 
disturbance near to the fast-surface. As a result of 
the deviation of waves, 1/Im(<5), and hence fluid am- 
plitude, rapidly decreases with 9 as indicated in Fig. || 
for a prolate shape of the fast-surface. We can quan- 
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equatorial plane 




Fig. 3. — Schematic figure (side view) of a black 
hole magnctosphcrc. The thick solid curves de- 
note the characteristics of equation (62). The fast- 
magnetosonic separatrix surface is denoted by the 
dashed line. The magnetic field lines are not drawn in 
order to avoid complication. The super-fast region is 
very thin in the magnetically dominated limit; how- 
ever, the width of this region is exaggerated. 



0.25 0.5 0.75 1 1.25 1.5 
polar angle 



Fig. 2. — Logarithmic variation of 1/Im(<5) as a func- 
tion of 9. 
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titatively understand this behavior by considering an 
approximate solution 



Im(<5) oc exp 



dxp 



(W 



(64) 



which is applicable when cPxp/dO 2 w 0. Equation 
(64) indicates that Im(<5) decreases exponentially as 9 
decreases. 

Let us finally consider the relation of amplitude 
among various perturbed quantities. Combining equa- 
tions presented in the last section, we obtain near to 
the fast-point 



Eg/M' B 



-,u 



£u r « u r 
hi 



(65) 



in order of magnitude because of k± m \JE/Jik r *. 
Other electromagnetic quantities have much smaller 
amplitude near to the fast-surface. It is interesting to 
note that the equipartition of energy ([e r ] 2 « /j,n[u r ] 2 ) 
is achieved near to the fast-surface, although the per- 
turbation energy is supplied mainly in the form of 
electromagnetic disturbances ([e r ] 2 « [E / /j]fxn[u r } 2 3> 
fj,n[u r ] 2 ) far from the horizon, In other words, a lot 
of perturbation energy is deposited from the electro- 
magnetic field to the fluid during the infall as a result 
of effective MHD interactions. In a realistic magnet- 
ically dominated black hole magnetosphere, the fluid 
becomes highly variable ([u r ] 2 ss [^ r ] 2 ) near to the 
fast-surface for a very small input of perturbation en- 
ergy ([e r ] 2 w [fi I E][Eg I M] 2 ) from the surroundings. 
The existence of the horizon is essential to make the 
fluid be highly variable near to the fast-surface (Pa- 
perl, and II). If the fast-surface is located far from the 
hole, fluid quantities does not become highly variable. 

So far, we have derived the following conclusions: 
(1) The fluid quantities become highly variable near 
to the fast-surface owing to MHD interactions near 
the horizon. (2) Their amplitude is a peaking func- 
tion at the separatrix surface located slightly inside 
of the fast-surface. This is because a large-amplitude 
fluid disturbance propagates along the separatrix as 
a outgoing fast-wave. (3) For prolate shape of the 
fast-surface, for instance, The fast-waves that prop- 
agate towards the equator can reach the fast-surface 
and escape to the sub-fast region. These results drive 
us to the question how these fast-waves propagate in 
the sub-fast regions. In the next section, we will be 
concerned with this issue. 



4. Wave propagation in sub-fast region 

4.1. Geometrical optics in relativistic accre- 
tion 

As we have seen, meridional wavelength is much 
less than the radius of curvature near to the fast sur- 
face. It follows that the propagation of wave packets 
of fast-magnetosonic mode follows the laws of geo- 
metric optics. In geometrical optics, the nature of the 
second-order partial differential equations describing 
the propagation can be well studied by the character- 
istic hypersurfaces of the system. The characteristic 
hypersurfaces, which play the role of wave fronts, can 
be expressed by a surface of ^{x^) = constant, where 
ip satisfies the following cikonal equation in the cold 
limit (Lichnerowicz 1967; see also Takahashi et al. 
1990 for sound waves): 



where s Q/3 is defined by 



EE 



U a UP 

TJ 2 ' 



(66) 



(67) 



Uym is the fast-magnetosonic speed in dr-basis and 
is defined by 



Ui 



FM — 



KB 2 p + B 2 /p w 2 



47r/in 
-KB r 



(68) 



The first term in equation (67) describes the influence 
of the gravitational field, while the second term that 
of cold, relativistic magnetohydrodynamic flows. If 
we were to replace ip i<x with k a , we would obtain the 
dispersion relation for the fast-magnetosonic mode, 



TT 2 h^h 



0. 



(69) 



Instead of solving the partial differential equation 
(|66|), we can investigate the trajectories of wave pack- 
ets by solving a set of the following ordinary differen- 
tial equations: 



dx a 
~d\ 
dp a 
dX 



dH(xP,pp) 

dH{xP, Pp ) 
dx a 



(70) 
(71) 
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where A is the parameter along a ray path. Since the 
Hamiltonian H contains neither t nor </>, both wave 
frequency lo — P t and azimuthal wave number m = 
—P<f> are conserved along a ray path. 

The unperturbed fluid's velocity field, on which the 
wave packets propagate, must be solved consistently 
with the equations of motion. First, the definition of 
proper time gives the poloidal wind equation 

„ (tttV i i _ ~ ZgyUtU* + 9tt(Ut) 2 

Pw 

(72) 

Secondly combining the unperturbed continuity equa- 
tion, Maxwell equations, and the frozen-in conditions 
with equations (^|) and (^), we obtain (Camenzind 
1986b) 



(gtt +g t 0^ F )e - M 2 E 

K-M 2 
{9t4, + ff00^F)e + M 2 L 



K-M 2 

where the Alfvenic Mach number M is defined as 



(73) 
(74) 



> 47T/X7? 2 



M 2 = 



B r 



(75) 



We assume an appropriate functional form for B r in- 
stead of solving the unperturbed trans-field equation. 

These three equations, ([72]) -([74]), together with 
(75), give the fluid's velocity field (U t , U r , —U^) 
on the poloidal plane. We assume that the accre- 
tion along each radial field line starts from the point 
at which K becomes 0.55. This condition defines a 
nearly spherical (but somewhat oblate) injection sur- 
face of accretion at n « 5M for a mildly rotating hole 
(a w 0.5M). We suppose that there is no flow of plas- 
mas outside of the injection surface. This assumption 
alters the propagation of the fast waves negligibly be- 
cause plasma flows in the region 5M < r < 10M is 
non-relativistic, whereas the fast-magnetosonic speed 
in di-basis is slightly smaller than that of light. We 
further assume that the energy density of the mag- 
netic field is nine times larger than that of the fluid's 



rest-mass in the sense that E/([xsm 9) 



TO. 



In this section, we trace the ray paths of the 
fast-magnetosonic wave packets radiated meridion- 
ally with momentum \k±\ ~ y ' E '/ ' \i |fc r *| from the 
fast-magnetosonic surface rather than radiated spher- 
ically (i.e., in all directions), by solving equations ( |70| ) 



and ( |7lD in the magnetically dominated accretion de- 
scribed by equations d7|)-(74). 



4.2. Collimation of MHD waves 

Let us first demonstrate typical results when the 
fast-magnetosonic speed, Ufm is slower in polar re- 
gions than in equatorial regions. Such a distribution 
of Ufm will be realized when \B r \ is smaller in the 
polar regions. A good example of such a magnetic 
field was presented by Blandford & Znajek (1977); 
they solved the vacuum Maxwell equations in a Kerr 
space time and derived a split monopole field, B r = 
_Bosin6> + 0(a 2 /M 2 ) for a distribution of a toroidal 
surface current density of I cx r~ 2 . Specifically, we as- 
sume that B r /(Awnr]) = {-E/n sin 2 0)(1 - 0.8cos<9) 2 
in this paper and calculate Ufm along a ray path by 
solving equation (pSS]). In this case, \rj\ becomes larger 
in the polar regions; therefore, the distribution of the 
fast-surface (Paper 1) 

- ( im f M0)w (76) 

r H (r H - MjMrn 

becomes prolate. Here, the function Ai(8) is of order 
of unity and has a weak dependence on 9 as 



At 



VKi-K 



(1 - aQ H sin 2 0)(1 - afl F sin 2 6) 3 W 3 



(77) 



where W is a function of 9 and of order of unity. For 
radial distribution of field lines, W 2 becomes 



1 + 



2ar H (^H - ^f) sin 2 9 



(r H - M)(l - afi H sin 2 0)(1 - aft F sin 2 9) 

(78) 

The information on the injection point of accretion 
appears only through K\ — 0.55. 

Examples of ray paths of axisymmetric waves (m = 
0) are presented in Fig. ||. All the wave packets are 
radiated from < 9 < ir/2 (the first quadrant) in this 
figure. Even though the wave packets have no angu- 
lar momenta (m = 0), they have non-zero angular 
velocities because of the space-time dragging and the 
rotational motion of the accretion flow. Therefore, 
the ray paths are projected on their instantaneous 
poloidal plane and are depicted in the figure. Since 
Ufm is smaller in higher latitudes, the fast-surface 
becomes prolate; thus the waves that propagate to- 
wards the equator can escape into the sub-fast re- 
gions. As a result, the wave packets radiated from 
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the first quadrant propagate clockwise as depicted 
in this figure. Because of the accretion, waves are 
pushed backwards to the hole and revolve around it. 
The heavy solid line on the equatorial plane denotes 
a dense disk which possibly resides around an active 
bole. If a wave packet collides with the disk, it will 
be totally absorbed to heat the plasma there. 

This figure indicates that most of the wave packets 
which are radiated meridionally from the fast-surface, 
collimate into the polar regions where J7fm * s small. 
Qualitatively the same results are obtained for non- 
axisymmetric waves. In Fig. [5], m = 2, 4, 6, 8, 10 waves 
are depicted. Ray paths in the first and the fourth 
quadrant originate at 9 = 30°, whereas those in the 
second and the third quadrant originate at 9 = —45° 
(in the second quadrant). Wave packets with negative 
angular momentum (m < 0) are not drawn, because 
they are soon swallowed by the hole. As the figure 
indicates, non-axisymmctric waves cannot reach the 
rotational axis because of their non-zero angular mo- 
menta; therefore, the tendency of collimation is some- 
what weakened compared with m = modes depicted 
in Fig. [|. 

In section 6, we discuss that the collimated waves 
may experience resonance and mode-convert itself 
into electromagnetic waves to result in a particle ac- 
celeration by non-linear interactions. Before we come 
to this issue, one more point must be clarified: if Ufm 
are larger in the polar regions than in the equatorial 
regions, ray paths must be bent towards the equato- 
rial plane. We examine briefly this case in the next 
subsection. 

4.3. Formation of a focal ring 

We demonstrate here typical results when the fast- 
magnetosonic speed, Ufm, is faster in the polar re- 
gions than in the equatorial regions. In this case, the 
distribution of the fast-surface becomes oblate. Ex- 
amples of ray paths for axisymmetric modes (m = 0) 
are presented in Fig. ||. All the wave packets are ra- 
diated from < 9 < tt/2 (the first quadrant) in this 
figure. Since the fast-surface is oblate, the wave pack- 
ets that propagate initially towards the rotational axis 
can escape into the sub-fast regions. Therefore, the 
waves propagate counterclockwise. 

This figure indicates that most of the wave pack- 
ets bent towards the equatorial disk to form a focal 
ring of radius ~ 5M and do not collimate towards 
the rotational axis. For non-axisymmetric waves, this 




/ 

/ 

/ 



'-10 -5 5 10 

radius (GM/c°2) 

Fig. 4. — Side view of ray paths of fast-magnetosonic 
wave packets in a magnetically dominated accretion 
flow around a mildly rotating black hole of a = 0.5M. 
Ray paths are radiated meridionally from the fast- 
surface every 5 degrees in the first quadrant: Solid 
curves denote ray paths radiated from the high lat- 
itudes (5°, 10°, 15°, 20°, 25°, and 30°). Dashed 
curves from the middle latitudes (35°, 40°, 45°, 50°, 
55°, and 60° degrees), while dotted curves from the 
low latitudes (65°, 70°, 75°, 80°, and 85°). Fast- 
magnetosonic speed is lower in the polar regions than 
in the equatorial regions. As a result, fast waves col- 
limate towards the rotation axis. 
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tendency is strengthened because of their non-zero 
angular momenta. 



Discussion 
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Fig. 5. — Side view of the ray paths of non- 
axisymmctric mode. The solid lines denote m = 2 
mode, while the dashed, dash-dotted, dotted, and 
dash-dot-dot-dotted lines denote m = 4, m = 6, 
m = 8 m = 10 modes, respectively. Waves propagat- 
ing in the right hemisphere are radiated meridionally 
at = 30° , whereas those in the left one are radiated 
at 6 = -45°. 



We have demonstrated that the magnetically dom- 
inated plasma accretion becomes highly variable near 
to the fast surface located close to the horizon and 
that such fluid's disturbance propagates as a fast- 
magnetosonic wave and collimates towards the ro- 
tational axis (especially for an axisymmetric mode) 
when the fast-magnetosonic speed is slower in the 
polar regions than in the equatorial regions. In the 
framework of MHD, the collimated fast waves will not 
cause interesting phenomena such as particle acceler- 
ation, even in nonlinear regimes. However, if we take 
the effects of plasma oscillation and cyclotron motion 
of particles into account, interesting results such as 
particle acceleration at a resonant point may be ob- 
tained (Holcomb and Tajima 1992; Daniel and Tajima 
1977). For this reason, we consider in this section a 
plasma wave in a pure electron-positron plasma (for 
observational and theoretical discussion on the exis- 
tence of electron-positron plasmas in AGN jets, see 
Ghisellini et al. 1992; Morrison et al. 1992; Xie et al. 
1995; Reynolds et al. 1996). 

As Fig. [| indicates, wave vectors (k) of colli- 
mated fast waves are nearly parallel to the magnetic 
held (B), of which poloidal components dominate the 
toroidal one near the rotational axis. It is, therefore, 
possible to regard that the collimated waves propa- 
gate as shear Alfven waves. 

In fact, a shear Alfven mode is one of the two pri- 
mary modes in a pure electron-positron plasma when 
k || B. The other mode is an electromagnetic mode. 
Their dispersion relation is given by 



to 2 (a 



2tV) 



c 2 (cj 2 - n c 2 ) 



(79) 



where the plasma and cyclotron frequencies are de- 
fined by 



p 



= ^4ime 2 /m e , (80) 
n c = — • (81) 

If we would set O e ^ u>, u> p , we could recover the 
dispersion relation of an usual Alfven mode, which is 
identical with that of a fast-magnetosonic mode when 
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k || B in a cold plasma. Exactly speaking, we must 
take account of geometrical correction due to hole's 
gravity in the definitions ( |S0| ) and (81); however, at 
the height where r > 10M, such effects become neg- 
ligible. 



From the dispersion relation ( 79 ) , it follows that a 
shear Alfven mode exists for frequencies less than O e 
(the resonance frequency), while an electromagnetic 
mode exists for frequencies greater than V^c 2 + 2cjp 
(the cut-off frequency) . In a realistic black hole mag- 
netosphere, f2 e is a decreasing functions of distance 
(r) from the hole, because B r decreases with r. 

It is important to note that the point of cut- 
off is located slightly outside of the point of reso- 
nance in a magnetically dominated magnetosphere 
(m e n e c 2 <C B 2 /8ir). Therefore, we can depict the fol- 



lowing scenario according to Budden (1961): A shear 
Alfven wave packet is injected outwards by the mecha- 
nism described in the preceding sections. As the wave 
approaches the point of resonance (a magnetic beach) , 
it increases the amplitude owing to a nonlinear effect. 
After reaching the point of resonance, it evanesces 
through the thin evanescent region to transmit as an 
electromagnetic wave above the point of cut-off. 

When the amplitude of injected shear Alfven wave 
is high and nonlinear effects become important, the 
wave may capture many particles to hide the singu- 
larity. In this case, they continue to propagate as 
an Alfven mode after passing the 'singularity' owing 
to the so-called self-induced tunneling (e.g., McCall, 
Hahn 1969). As a result, the trapped particles may 
be accelerated through long distances to become rel- 
ativistic. According to the nonlinear simulations per- 
formed by Daniel and Tajima (1977), who adopted 
a very thin evanescent layer corresponding to a mag- 
netic dominance of B 2 /8tt ~ 9n c m e c 2 , particle accel- 
eration up to energies of 8m e c 2 can be realized if the 
injected wave is highly nonlinear. On these grounds, 
it seems reasonable to conclude that the process of 
wave amplification and collimated demonstrated in 
this paper triggers initial acceleration of jet due to a 
cyclotron resonance in an electron-positron plasma. 
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Fig. 6. — The same figure as Fig. 4; however, the 
fast-magnetosonic speed is slower in the equatorial 
regions than in the polar regions. As a result, fast 
waves are bent towards the equatorial plane to form 
a focal ring. 



Thanks are due to Drs. A. Tomimatsu, K. Shibata, 
and M. Takahashi for valuable advice and helpful sug- 
gestions. 
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